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Abstract 

In the spirit of previous papers, but using more general field configurations, 
the non-linear 0(3) model in (2+l)-D, modified by the addition of both a 
potential-like term and a Skyrme-like term, is considered. The instanton so- 
lutions are numerically evolved in time and some of their stability properties 
__|- studied. They are found to be stable, and a repulsive force is seen to exist 



among them. These results, which are restricted to the case of zero-speed 

in 
23 



systems, confirm those obtained in previous investigations, in which a similar 
problem was studied for a different choice of the potential-like term. 



%*' 1 Introduction 

In the past few years, a-models in low dimensions have become an increasingly im- 
portant area of research, often arising as approximate models in the contexts of both 
particle and solid state physics. They have been used in the construction of high-T c 
superconductivity and the quantum Hall effect; in two Euclidean dimensions, they 
appear to be the low- dimensional analogues of four-dimensional Yang-Mills theories. 
But only very special cr-models in (2+l)-D are integrable and the physically rel- 
evant Lorentz-invariant models are not amongst them; in these cases, recourse to 
numerical evolution must be made. The simplest Lorentz-invariant model in (2+l)-D 
is the 0(3) model, which involves three real scalar fields, (^(x At )={c J > a (x A1 ), a=l,2,3}, 
with the constraint that <f) lies on the unit sphere S2 ■ 

0.0=1. (1) 

Subject to this constraint, the Lagrangian density and the corresponding equations 
of motion are 

£ a = \(d,cf>).(d»cf>), (2) 
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d 2 4> = [-(d t <p) 2 + (<9,</>) 2 + (d y( j>) 2 }cj> + dl4> + d 2 y <f>. (3) 

Note that we are concerned with the model in (2+l)-D: x^ = (a; , x 1 , x 2 ) = (t,x,y), 
with the speed of light set equal to unity. 

An alternative and convenient formulation of the model is in terms of one inde- 
pendent complex field, W, related to the fields a via 

w = ±rir- « 

01 + 202 

In this formulation, the Lagrangian density and the corresponding equations of motion 
read (asterisk denotes complex conjugation) 

d,Wd»W* 

(1 + |W| 2 ) 2 ' [b) 



and 

2W*[(d t W) 2 - (d x W ) 2 - (d y W) 
l + \W 



d 2 w = d 2 x w+d 2 w+ ^ y [y t ' 7:;z: ( 6 ) 



The problem is completely specified by giving the boundary conditions. As usual we 
take 

lim cl>{r,6,t) = <j> (t), (7) 

where <f>o(t) is independent of the polar angle 9. In (2+l)-D this condition ensures 
a finite potential energy, whereas in two Euclidean dimensions, i.e., when </> is in- 
dependent of time, this leads to the finiteness of the action, which is precisely the 
requirement for quantization in terms of path integrals. As shown by several people 
0. ||, any rational function W(z) or W(z*), where z = x + iy, is a static solution 
of eq. (^). These are the instantons of the model, and can be regarded as static 
solitons of the same model in (2+l)-D. The simplest one-soliton solution, W = Xz 
(A is a free parameter determining the size of W), has been numerically studied by 
Leese et.al [|J. When viewed as an evolving structure in (2+l)-D, the soliton has 
been found to be unstable; any small perturbation, either explicit or introduced by 
the discretization procedure, changes its size. This instability is associated with the 
conformal invariance of the 0(3) Lagrangian in two dimensions. 

The 0(3) solitons, however, can be stabilized through a judicious introduction 
of a scale into the model, thereby breaking its conformal invariance. As explained 
below, in the present article we study this problem following the methods of a pre- 
vious investigation ||, but considering a different potential-like term and restricting 
ourselves to the systems initially at rest. 
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2 Skyrme model in (2+1) dimensions 



Using the I^-formulation, our Skyrme model is defined by the Lorentz-invariant La- 
gragian density 

C = C a 

nQ Ad t W*d y W - d t Wd y W*) 2 + (d t W*d x W - d t Wd x W*f 

~ i[ {i+ww 

{d x w*d y w -d x wd y w*f 
{i + \w\ 2 y J 

\w 

- A9o 





w 


-A| 8 


(1+1 


W\ 


2)4 



where L a is given by eq. (^) and #i,#2 are real parameters with dimensions of length 
squared and inverse length squared, respectively; they introduce a scale into the 
model, which is no longer conformal invariant. If the size of the solitons is appropri- 
ately chosen, it is energetically unfavourable for the solitons to change it. The ^i-term 
is the (2+l)-D analogue of the Skyrme term, whereas the #2-term is a potential-like 
one. Unlike the former, the latter term is highly nonunique [[J. 

The field equation corresponding to the above Lagrangian can be cast into the 
form 



w tt = w xx + w vv + 



2W*[(W t ) 2 - (W x ) 2 - (W t 



\2l 



yy ' 1 + \W\ 2 

A6 1 



1+ W 



—[2w; x w t w x + 2\v; y w t w y - 2\v* v w x w y 



+ w* xx (w 2 - w 2 ) + w; y (w 2 - w 2 ) - w; t (w 2 + w 2 ) 

+ W XX (\W t \ 2 - \Wy\ 2 ) + Wyy(\W t \ 2 ~ \W X \ 2 ) + W tt (\W X \ 2 + \Wy\ 2 ) 

+ w X y{w* x w y + w;w x ) - w tx (w;w x + w* x w t ) - w ty (w;w y + w;w t ) 



l + \W\ 
169 2 \W - X 
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+ (l + \W\ 2 )* ' (9) 

where the notation W x = d x W, W xx = d^W, etc., has been used. 
It is straightforward to check that 

w = \ z -^4 (10) 

z — b 

is a static solution of eq. (0), provided the following relation holds: 



A = ^r^- (11) 
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Now A, which characterizes the size of the instanton W, is no longer a free pa- 
rameter. It is fixed by eq. The instanton, with its size thus fixed, is usually 
referred to as a 'skyrmion'. Unlike the ordinary 0(3) case [|J, where the instanton 
changes its size as time elapses, we expect our skyrmion to be stable. In reference 
5, a similar question was studied using —402/(1 + \W\ 2 ) 4 for the potential-like term, 
with the Skyrme field given by Xz, where A = d26\/ 02- 

3 Numerical setup and results 

To perform our numerical simulations we employ the fouth-order Runge-Kutta 
method, and approximate the spatial derivatives by finite differences; the Lapla- 
cian is evaluated using the nine-point formula. All computations were performed on 
the workstations at Durham, on a fixed 201x201 lattice with spatial and time steps 
5x=5y=0.02 and 5t=0.005. Every few iterations we rescaled the fields — > <f>/ \J4>-4> ■ 
We also included along the boundary a narrow strip to absorb the various radia- 
tion waves, thus reducing their effects on the skyrmions via the reflections from the 
boundary. As time elapses, this absorption of radiation manifests itself through a 
small decrease of the total energy, which gradually stabilizes as the radiation waves 
are gradually absorbed. 

We choose the parameters to have the values: #i=0. 015006250, 6> 2 =0.1250, a=0.75 
and 6=0.05 which, according to eq. (0), set A=l. 

3.1 Results for one skyrmion 

The numerical evolution shows that the skyrmion's shape, of which a representative 
picture is given in the upper half of figure 1, practically does not change as a function 
of time. The maximum of the total energy density, which is related to the soliton size, 
starts off at the value 128.5 and, after some radiation waves are emitted, stabilizes at 
~ 129.5 (see lower half of figure 1). This is in agreement with the analytical result, as 
the expression for the maximum of the total energy density, namely, E max = e(l+0\e), 
£ = 8(| A | 2 + 1) 2 /| A(a - b) | 2 , yields 129.3, the 'canonical size'. 

In figure 2 are exhibited some kinetic energy density waves. They propagate out 
to the boundary at the speed of light, leaving the region where the skyrmion is located 
essentially free of kinetic energy. The smallness of the kinetic energy indicates that 
the field given by eq. ( |i~0"D remains almost perfectly static. For example, the amplitude 
of the total energy density, at t=0, is positioned at z max =(0, 0.4000) and, for t=10, 
it has been slowly shifted to (0, 0.4013). 

So, as expected, our skyrmion appears to be stable under perturbations brought 
about by the discretization procedure. These results agree with those of reference ||, 
where the skyrmion was positioned at the centre of the grid, i.e., z max =0, which, as 
one can easily check, corresponds to a skyrmion given by W = Xz or W = 1/Xz. 
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In order to study the stability problem further, we are currently performing simula- 
tions with the same ^-values but setting A ^ 1. This corresponds to initial conditions 
that are no longer those of the static solution. We hope some concrete results will be 
available in the future. 



3.2 Results for two skyrmions 

Let us now consider the case where two instantons are put on a lattice. We consider 
the field 

w = x — h ^ < 12 > 

z — o z + a 

with A=l, a = c = 0.75, b = d = 0.05 and let it evolve. One can readily verify that 
the above field configuration is not a solution of the equations of motion, and hence 
it evolves when started off at rest. 

Our simulations have shown that, at the initial time, E max is well above the canon- 
ical size as determined in the last sub-section. As soon as the evolution in time begins 
the solitons shake off some kinetic energy and alter their size by getting broader; E max 
decreases and then undergoes damped-oscillations around a value very close to the 
canonical one. As this happens, the solitons slowly move away from each other, ex- 
hibiting the presence of a repulsive force between them. By t ~ 8, the amplitude of 
the above-mentioned oscillations is quite small, and E max eventually stabilizes near 
the expected canonical value. Typical pictures for this case are displayed in figures 3 
and 4. 
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Figure 1: Total energy density for one skyrmion at the initial time and the evolution 
of its peak. 



Kinetic Energy Density 

Kmax=0.000132; t-0.5 




Kmax=0.0000438; t=1.5 




Figure 2: Kinetic waves travelling towards the boundary. 
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Total energy density 

Emax=268.4; t=0.0 




Figure 3: Two skyrmions. A repulsive force exists between them. 
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Time 



Figure 4: Evolution of the maximum of the total energy density for the two-skyrmions 
case. 
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4 Conclusions 



We have performed a numerical study of the time evolution of the general one- 
instanton solutions of our version of the Skyrme model in (2+l)-D, confirming some 
of the properties previously seen in other versions of the model. The extra terms 
added to the ordinary non-linear 0(3) model break its conformal invariance and sta- 
bilize the solitons. The changes of soliton size and position are negligible as time 
progresses. Also, the time evolution of the two-skyrmion configuration has showed 
us that, like in the previously studied model, the forces between the skyrmions are 
repulsive. 

A paper on the scattering properties of our skyrmions will be available in not too 
distant a future. 
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